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C " 3 ' Abstract. We present doubler-free gauge-invariant lattice vector gauge ac- 

C " 3 , tion for some real representations of Wilson gauge fields on an octet of fermions. 

^Nj 1 It is based on a geometric representation of the Dirac equation as an evolution 

equation on the three-dimensional exterior bundle A(1R'^) for a single bispinor 
and of the bundle (A(g)A)(lR^) for an octet. We find doubler-free lattice Dirac 
) , operators for above bundles. A gauge-invariant connection with Wilson lattice 

■ gauge fields is possible for some real representations of the gauge group. The 

QCD action of 311(3) is of this type. 

Application in lattice QCD seems useful: We don't have to waste time and 
memory for doublers as well as for correction terms to suppress them. 
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1. Introduction 



One of the main problems of lattice gauge theory is the fermion doubling prob- 
lem. If we use a naive central differences discretization of the Dirac equation, we 
observe a doubling effect: The large distance limit of the resulting lattice equation 
i does not describe the original single Dirac particle, but, instead, sixteen Dirac par- 

ticles — a factor two doubling in each spacetime dimension. Such a discretization 
Oh, is inappropriate for quantum computations where the number of degrees of freedom 

is important. 

A partial solution is given by a geometric variant of the Dirac equation — 
the Dirac-Kahler equation [4]. This equation is defined on the complexified four- 
^> \ dimensional exterior bundle C (8> A(M'*). This bundle contains sixteen complex 

^ • components, and the Dirac-Kahler equation on it is equivalent to the classical 

5^ I Dirac equation for four Dirac bispinors. On the other hand, for this equation exists 

a doubler-free so-called "staggered" or Kogut-Suskind discretization ([5], see also 
[6]). Instead of sixteen doublers of the naive discretization this gives a lattice equa- 
tion for four Dirac bispinors. Another solution is to suppress the doublers using 
some additional terms which give all the doublers except one a large mass. The 
most famous one is the Wilson mass term [7], but many different versions of such 
mass terms have been proposed. 

In this paper we consider another variant of geometric fermions. Instead of 
the bundle C (8) A(R^) we use the three-dimensional exterior bundle A(R^) and a 
direction-dependent complex structure on it. The four-dimensional Dirac equation 
appears as an evolution equation on this bundle. The time t is only a parameter 
of dynamics without geometric meaning. This representation of the Dirac equation 
describes already a single bispinor. 

The main advantage of the geometric Dirac-Kahler fermions remains — there 
exists a doubler-free "staggered" discretization of the continuous Dirac equation. 
But now it gives a doubler-free "staggered" discretization of a single bispinor. We 
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find also a remarkable, more symmetric discretization for an octet of bispinors 
related with the bundle (A (g) A)(M^). 

Using the standard Wilson approach to lattice gauge fields becomes problem- 
atic because the complex structure does not act pointwise. But to put the octet 
representation of SU{3) used in QCD on the lattice we do not need the complex 
structure: It is equivalent to a real representation. As a consequence we obtain a 
doubler-free discretization for lattice QCD. 

The discretization presented here is closely relatcid with the "spacc-geomctric" 
interpretation of standard model fermions as sections of the bundle (T® A(g)A)(IR^), 
or, on curved background with ADM decomposition = © M, the bundle 
(T A O A)(M3) made in [8]. 



2. Space-geometric representation of spinors 

The classical Dirac algebra is defined by the operators 7* with commutation 
relations {7', I-'} = rj'^-' — the Clifford algebra CZ(1,3). But there are also some 
other operators which play an important role in the standard model: The complex 
structure, which is defined by the operator i, and a real structure defined by the 
operator C, which is necessary to define Majorana mass terms. The commutation 
relations of this "extended Dirac algebra" {7'', i, C} are: {7'*, 7^"} = 2r]^'^, [i, -y^] = 
0, [C, 7^^] = 0, {C, i} = 0, C2 = 1, = -1. 

The three-dimensional or space-geometric representation of spinors based on the 
following isomorphism: 

Theorem 1. The "extended Dirac algebra" {'y'^,i,C} is isomorph to the Clifford 
algebra CZ(3,3). The Clifford algebra generators are the 7' for i > and the three 
operators defined by 



(3' = C, 
/3' = Ci, 

^3 ^ ^5 ^ _ .^oy^2^3_ 

Proof: An explicit check proves that the operators 7*,/3* fulfil the properties 
required for the generators of Cl{3,3) (they anticommute, = —(7*)^ = !)• 

We have already expressed the generators of CZ(3, 3) in terms of the extended 
Dirac algebra {7'', i,C}. The reverse expressions can be easily found: 



(1) c = 

(2) i = 

(3) 7° = PWP'l'P'l' 

It remains to prove that there are no other relations in {7'*, i, C}. This can be 

done by an explicit computation of the dimension of the extended Dirac algebra. We 
have six generators {7^, i,C}, for each pair of generators we have one commutation 
resp. anticommutation rule, the square of all generators is ±1, and there are no 
other independent relations. Therefore an arbitrary monom may be transformed 
into a monom where every generator appears at most once. There are 2^ such 
monoms, which gives the dimension of C/(3,3). qed. 
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The Clifford algebra Cl{3, 3) acts in a natural way on the three-dimensional Clif- 
ford bundle CZ(M^). This bundle is equivalent to the exterior bundle A(R^) which 
consists of inhomogeneous differential forms ip = (fin{x){dx^)'^^{dx'^)'^^{dx'^)'^^, 
K = (ki,K2,K3) G (^2)^. The operators 7* act as in standard Hodge theory and 
define the Hodge theory Dirac operator 

(4) = ij% 

(see, for example, [1]). The action of the operators /3* follows from the following 
formula: 

(5) (/3Sv)«w = i-irM^) 

It is not hard to see the following: 

Theorem 2. The representation of the extended Dirac algebra {7'', i, C} on A(IR^) 
defined by (4), (5) and theorem 1 is equivalent to the standard representation. 

Theorem 1 gives now the following four- dimensional Dirac operator: 

(6) Di = n"9o + 

The definition of the Lagrange formalism simplifies if we use the following variant 
of the Dirac operator: 

(7) D = -i7S%^£'4 = idt + ta'di 

where i = /3^/3^/3^, = —1, {<•, 7'} = 0, [i, a'] = 0. The equivalence is obvious: 
Dip = iff D4(p = 0. Now wc can use the standard Euclidean scalar product on 
A(M'^) to define the Langrange formalism as 

(8) L = {ipD^) 

The Dirac operator considered here should not be mingled with the four- 

dimcnsional Dirac-Kahlcr operator defined on the spacctimc bimdlc C ® A(R'*) 
which has real dimension 32. It defines an evolution equation on the spatial bundle 
A(M^) with real dimension 8. Time is only a parameter of dynamics and has no 
geometric interpretation. Manifest Lorcntz invariancc is lost. Only the equivalence 
theorem proves that Lorentz invariance will be recovered for observable effects. 

This representation also violates manifest symmetry in space: The complex 
structure i = P"^ 0^ as well as the operator = fi^ depend on a preferred di- 
rection in space. The introduction of three generations of fermions associated with 
the three spatial directions may be used to get rid of this asymmetry. 

A similar " space Clifford" representation of the Dirac equation has been consid- 
ered by Daviau [2], see also [3]. 

A comparison of space-geometric spinors with standard spinors can be found 
in [8]. It is shown there that on Minkowski background we have no observable 
differences, but on curved background we obtain differences which are observable 
in a graviational analogon of the Bohm-Aharonov experiment. 
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3. Lattice theory for the three-dimensional Dirac operator 

Let's consider at first the naive central differences discretization of the Dirac 
operator D on A(R^). We use a regular rectangular lattice with lattice functions 

€ ri : Z| (g) Z'' ^ R denoted by ipi^{n) € R, k = {ki,K2,k^) G (^2)^, n = 
(7x0,711,712,^3) G Z"'. Using the denotations = {5o^,5i^, 82^^53^) we define 



(9) {t''^'p)n{n) = ip^{n±h^) 

(10) = lif^^-tt) 

(11) Dh = iAo + La'Ai 

The Lagrange formulation of the lattice theory can be given in the same form 
as (8), but with the lattice version (11) of the Dirac operator D: 



(12) L = {^Dh^) 

3.1. Decomposition into sixteen staggered fermions. A well-known problem 
of lattice gauge theory — the so-called "fermion doubling problem" — is that 
the naive central differences discretization of the classical Dirac equation does not 
describe a single Dirac fermion, but sixteen. This doubling problem appears in our 
naive discretization (11) too. As usual, we obtain sixteen Dirac fermions instead of 
one. But in our discretization these sixteen doublers decompose in a quite simple 
way. 

Theorem 3. There exists a decomposition 

A 

A = (Ao,Ai, A2, A3) e Z^, Aj = ±1 so that for = Y.^p\, fx e O-^ 

DhV = ^ Dh^x = 

and the subspaces fl^ cSl are realized as staggered sublattices. 

To see this let's consider the following lattice operators: 



(13) {eo^).{n) = (-l)"''+"^+"^+">«(n) 

(14) (e.^).(n) = (-l)'='+"'^«(7x) 

These four operators commute with each other and fulfil = 1. This allows 
to define common eigenspaces with eigenvalues ±1. Thus, we define the fl^ as 
eigenspaccs e^Lp = \^(p. It is easy to verify that the Q,^ define staggered sublattices: 
They contain only the component defined by (—1 = Aj on nodes which fulfil 

the property (-l)"o+ni+n2+n3 = Xq. 

Now, the also anticommute with Dh. This gives the property : fl^ — » Q,~^. 
As a consequence, from Dhf = follows D^fx = for all A. 

The decomposition on the level of lattice components suggests to use one of the 
eigenspaces, say fif^'^'i'i), as a staggered discretization of a single Dirac fermion. 

A comparison of such a single staggered fermion with the standard staggered 
Dirac-Kahler fermions seems interesting: In the classical staggered fermion concept 
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[5] there remain four doublers (more accurate, already the continuous Dirac-Kahler 
equation describes four Dirac bispinors, it's staggered discretization is doubler- 
free as a discretization of the Dirac-Kahler equation but nonetheless describes four 
bispinors). Instead, in our approach we have no doublers, but a single Dirac fermion 
(eight real components). On the other hand, we have to pay for this with two 
disadvantages: We have no longer manifest, explicit relativistic symmetry, but a 
special role of the time coordinate. On curved background, such a special role for 
a time coordinate requires an additional structure an ADM decomposition. 

The other problem is the complex structure. First, it depends on a direction in 
space. This dependence disappears if we introduce throe generations of fermions 
and associate the three gcncirations with directions in space (see [8]). A more seri- 
ous problem is that the complex structure on the lattice is not pointwise: real and 
imaginary parts of the four complex components of a bispinor are located on differ- 
ent lattice nodes. Especially there is no pointwise operator for multiplication with 
i on our staggered fermion. Now, multiplication with i is obviously an important 
operation in the standard model. Even if we get rid of it in the pure Dirac equation 
using the form (6) or (7), they reappear in the interaction terms with gauge fields 
as well as in the mass terms. Fortunately for mass terms it seems sufficient to fix 
some non-pointwise replacement for i. For gauge fields, the non-pointwise character 
of multiplication with i seems more problematic. 

3.2. The geometric octet. To define the standard model (inclusive right-handed 
neutrinos), we have to consider three octets of Dirac fermions. Now, a straight- 
forward way to do define an octet would be to use eight copies of the staggered 
fermion defined by fl^^'^'^'^h But there is a much simpler and more natural way to 
obtain an octet of fermions which we name geometric octet: We use the eigenspace 
n C defined by eQ(p = ip which consists of all components on "even" nodes. The 
lattice space has not only a much simpler definition (we don't even have to de- 
fine the operators ej,i > to define Cl), but also a larger symmetry group: While 
on different lattice nodes in ri^^'^'^'i) we have different components all lattice 
nodes in f2 are of the same type: They contain all eight components 

Thus, from point of view of symmetry between lattice nodes, the geometric octet 
of Dirac fermions seems much more natural than a single Dirac fermion. Let's con- 
sider its geometric nature: The eight doublers ri(i>*J''=) in this octet are connected 
with each other by translations in space. The components of the continuous limit 
of the geometric octet therefore should be connected by infinitesimal translations. 
This gives transformational properties as for the geometric bundle (A ® A)(M^) 
instead of A(M3)8.i 

There is one problem with the reduction f7 ^ f7 - the naive Lagrange formula- 
tion (12) is not well-defined on Ct. But there is a simple solution — instead of (12) 
we can use the following Lagrangian: 

(15) L = i^t^Dhif) 



Of course, as long as we consider only the Minkowski background, as we do in this paper, 
there is no physical difference between (A® A)(K'^) and A(K'^)**. The difference becomes important 
on curved background. The consideration of discretizations for curved background we leave to 
future research. 
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where the translation defines a preferred direction in space. This requires the 
same type of symmetry breaking as the choice of the complex structure i = (3^/3^, 
and we can get rid of it if we introduce three fermion generations (see [8]). 

Below we focus our interest on the "geometric" octet. 

4. Lattice gauge theory 

Our lattice Dirac equation for the geometric octet may be connected with gauge 
fields following the classical Wilson scheme for lattice gauge fields, even with some 
simplifications. In the Wilson approach, the lattice gauge field for a gauge group 
G is described by gauge transformations Ufi,{n) G G defined on edges (links) of 
the lattice. The lattice gauge transformations are defined by functions g{n) G G 
on the nodes with values in the gauge group G. In comparison with standard 
lattice fcrmions we have no pointwisc operation of multiplication with i. Hut if the 
representation T : G ^ gl{S) which acts on the space spanned by the ^'^(n) is real, 
we can nonetheless use the standard way to define a gauge-invariant interaction 
Lagrangian — to replace the difference operators in the Lagrangian by covariant 
difference operators A^: 

(16) (A^^),(n) = ^(T{U^(n))^An + h^) - T(U-\n - h^.ji))^,{n - h^)) 

The restriction to pointwise gauge transformations (and the resulting restriction 
to real representations) is is not only prescribed by the standard Wilson approach. 
It also allows to preserve the reduction ^ O and, therefore, allows to define the 
gauge action not only on the original sixteen doublers but also on the geometric 
octet Cl. This gives the Lagrangian 

(17) L = {<fit^Dhip) + ■'-'Wilson 
for the gauged Dirac operator Dh = (-Aq + i-a'Aj. 

4.1. QCD on the geometric octet. Let's consider now the specific properties of 

local gauge actions on the geometric octet. To define QCD on the geometric octet 
we need the appropriate octet representation of the Lie algebra sm(3) in terms of 
pointwise operators which commute with the Dirac operator Dh- For this purpose, 
let's consider the following operator algebra: 

(18) A = 

(19) 7i = Ae^ 

We have /?? = 1, 7? = —1, and all Pi,ji anticommute, thus, it is the Clifford 
algebra CZ(3, 3) = gl{8). The compact subalgebra 0(8) C gl{8) is generated by 
products g of the generators which fulfil g"^ = —1. 

Moreover, we have simple commutation relations with the Dirac operator: 
[PijDh] = 0, {7j,Z)/t} = 0. If we require [g,Dh] = 0, we have to restrict our- 
self to monoms with even number of factors 7^. This restriction is equivalent to 
[g,i] = with i = P1P2P3, = —1- The subalgebra of 0(8) which fulfills this 
restriction is equivalent to w(4), where i plays the role of the related complex unit. 
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It should be noted that this is a real eight-dimensional representation of u{4). A 
complex unit i which acts on single fermions ri(i'*J''^) has not been used here. Such 
a complex unit i should be distinguished from the real operator i G gl{8) which 
does not preserve the decomposition into fermions fi^^'*'-?''^). 

As a real representation, this octet representation of u(4) would be irreducible. 
But as a complex octet representation it is reducible: We have an involution op- 
erator a, {iiY = 1, [g,ii] = for g e u(4). This allows to define two invariant 
subspaces with projectors = |(1 ± ii). On these four-dimensional subspaces 
fi± = P^ri wc have i = ^i? The observable quarks and leptons will be, therefore, 
linear combinations of pairs of staggered fermions and defined 

by the projectors P^. It seems natural to identify the pairs Q(i'»>J'''=) qqCI'-'.-J'-^) 
with clcctrowcak doublets. 

Thus, we obtain two times the defining four-dimensional representation of u(4). 
Restricted to su(3) C u(4) this gives two times the defining representation of su(3) 
and two times a singlet representation of sw(3), as required for the QCD action. 
At least as long as we do not consider weak interaction the two subspaces f2± are 
independent. Therefore the different choice of the sign of the complex structure is 
a question of convention. 

It remains to define an appropriate restriction from u(4) to sm(3). For this 
purpose we have the interesting symmetric operator 

(20) C7 = /?i7i/3272 + /?272/3373 + /3373/3i7i 

which defines baryon resp. lepton charge in a natural way. The restriction 
[g, a] = defines the subalgebra © w(3) C u{4). Its generators are: 

(21) 1 = 010203, ^i = 0il3lk, i^i= PjPk-^j^k, Vi=7i{Pj^j - Pk7k), 

with l<i<3, j = i + l mod 3, A; = i -|- 2 mod 3. The two-dimensional center 
of this algebra is generated by the generators i, S = 5i + 82 + S3. We obtain su{3) 
as the remaining factor algebra, with the eight generators oji, rji, Si — S2, S2 — S^. 

4.2. The mass terms. To put various mass terms on the lattice we need lattice 
versions of the complex unit i = f3^P^ as well as of the operators C = /3^, 7^ = 
These have to be operators which leave the siibspaces fi'*' invariant. Because they 
are transformations of different components (f^, which are located on different nodes 
in the staggered subspaces fl^, such operators cannot be pointwise operators. As 
a consequence, it is not possible to preserve all important exact properties of the 
continuous operators i,C, 7^ on the lattice. The question is which of the exact 
properties are the more important. For example, a variant for i which commutes 
with Dh is: 

(22) iih^)^in) = (/3VV)«(n + hi + /12). 

But this operator does not define a complex structure, (ih)^ 7^ —1 but includes 
a shift. We can, instead, preserve the property = —1, for example with: 

■^Notc that on the lattice such a decomposition would require a lattice aiialogon of i with exact 
properties = —1, = 0, [i,-D/i] = which does not exist. Thus, this decomposition will not 
be exact on the lattice. It appears only in the continuous limit. 
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(23) {ih^).{n) = + (-l)"^/i2) 

In this case [ih,Dh] ^ 0. Similar variants exist for the operators C, 7^. Which 
choices for these operators are optimal for the implementation of mass terms has 
to be left to future research. 

5. Discussion 

Wc have defined a doublcr-free gauge-invariant lattice gauge theory for the SU{3) 
octet action of QCD. Such a discretization is in itself quite useful for applications in 
lattice QCD: In comparison with the standard approach to lattice QCD we do not 
have to waste time and memory for the doublers, and we do not have to compute 
the additional terms to suppress them. 

But this discretization seems interesting not only because of these applications. 
It has been found by consideration of a "space-geometric" interpretation of the 
whole fermionic sector of the standard model as the bundle (T(g)A(g)A)(IR^) proposed 
in [8]. The vector index describes the three fermion generations, and the bundle 
(A (E) A)(R''') is the largo distance limit of our geometric octet. The direction- 
dependent structures of our octet (the complex structure i = /3^/3^, use of in 
(15)) may be combined into direction- independent structures if we consider all 
three generations and associate these generations with different preferred directions 
in space. 

On the other hand, there is no simple way to recover manifest relativistic symme- 
try. In the; space-geonicitric interpretation time is only a non-geometric parameter 
of evolution. Geometric interpretation is only three-dimensional. This is in conflict 
with the four-dimensional spacetime interpretation of the gravitational field, but in 
natural correspondence with the ADM decomposition a foliation of spacetime 
into space and time = © R. In an ADM decomposition the gravitational 
field is described by a positive-definite three- metric, a three-dimensional vector 
field, and a positive scalar field. In [8] we propose to interpret standard model 
fermions as sections of the bundle (T (g) A (g) A)(M^). 

If this discretization really has some fundamental importance, then there should 
be also a natural possibility to describe electroweak interaction on the geometric 
octet. The key for understanding chiral gauge fields may be a better understanding 
of the necessarily non-pointwise lattice approximations of i, C, and 7^. Note, for 
example, that the use of an ih which preserves {ih)^ = —1 like (23) defines a de- 
composition into two lattice subspaces which is not exactly preserved by Dh- A 
decomposition into exact mass eigenstates on the lattice will give different subspaces 
which will not be left invariant by the QCD gauge fields. They will coincide only 
in the continuous limit. What will be the highest order correction terms? It will 
be some interaction between the mass eigenstates inside the electroweak doublets 

Q(l,i,j,fe) Q ^(l,-i,-j,-fe)_ 

Maybe for some appropriate choice of the mass terms this discretization effect 
gives electroweak interaction or something sufficiently close to it? This hypothesis 
we leave to future research. 
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